Collapse Dynamics of a Polymer Chain: Theory and Simulation 
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We present a scaling theory describing the collapse of a homopolymer chain in poor solvent. At 
time t after the beginning of the collapse, the original Gaussian chain of length N is streamlined to 
form N /g segments of length R{t), each containing g ^ t monomers. These segments are statistical 
quantities representing cylinders of length R ~ t^^'^ and diameter d ^ t^^"^, but structured out 
of stretched arrays of spherical globules. This prescription incorporates the capillary instability. 
We compare the time-dependent structure factor derived for our theory with that obtained from 
ultra-large-scale molecular dynamics simulation with explicit solvent. This is the first time such 
a detailed comparison of theoretical and simulation predictions of collapsing chain structure has 
been attempted. The favorable agreement between the theoretical and computed structure factors 
supports the picture of the coarse-graining process during polymer collapse. 

PACS numbers: 36.20.-r,64.60.Ak,64.60.Fr,61.25.Hq,83.10.Nn 



The collapse transition of a single polymer molecule 
is for several reasons an interesting physical problem.]^, 
|, |, |, |, |, 0, |, g Understanding the collapse tran- 
sition is foremost a precursor to understanding pro- 
tein folding, px|] and how other biomolecules react to 
changes in environment. Furthermore, a better phys- 
ical picture of the collapse transition will enhance the 
ability of nanoscale micromanipulation techniques to 
provide a way to experimentally approach many such 
problems.!^, |ll| 

The earliest consistent theoretical picture of the col- 
lapse transition stems from de Gennes.||l| In this sce- 
nario, when the temperature is shifted by AT from 6 
conditions, a Gaussian coil starts to aggregate, resulting 
in a uniformly dense sausage-like shape. As time goes on, 
the minimization of interfacial area drives this sausage- 
like shape to thicken and shorten until, at the final stage, 
a globule is formed. More recent work has shown that 
such uniform sausage-like shapes are highly unstable in 
solvents due to the capillary instability, which selects for 
"pearl- necklace" structures over uniform sausages.]^ At- 
tempts made to incorporate the capillary instability in 
the description of the collapse of a polymer chain, [p2|, 
however, have so far resulted in no universally accepted 
theoretical description, particularly in prediction of how 
the time for collapse scales with chain length. The pur- 
pose of this paper is to present a theory based on a new 
idea for the mechanism of the collapse transition, and to 
support this theory by comparing directly to results of 
large-scale molecular dynamics simulations. 

Consider a Gaussian coil representing a polymer in 9 
solvent. We may think of this coil as a Gaussian fractal, 
created by a heirarchical scheme, demonstrated in the 
right side of Fig. ^ We subject this system to an instan- 
taneous drop in temperature, which quenches the system 
rapidly into poor solvent conditions. This effectively in- 
troduces monomer-monomer attraction. We assume this 
quench is so deep that the initial nucleation of the first 
globules of the collapsed polymer phase appear instan- 




FIG. 1: 

fractal. 



Equivalence between a Gaussian coil and a Gaussian 



taneously. (Otherwise, addition time scales and regimes 
would enter the problem.]^, |l^, The system evolves 
dynamically according to microscopic equations of mo- 
tion of chain segments. The segments are straightened 
arrays of globules, and further straightening streamlines 
them at larger length scales. This leads to the overall 
coarse-graining of the chain, on a length scale that is 
growing with time. This process is the reverse of the 
heirarchical Gaussian fractal generation. 

Because we view the original Gaussian chain made of N 
monomers of size , we can define a length scale R as the 
resolution at which the chain appears as a Gaussian chain 
made of N / g coarse-grained monomers, or segments, of 
size b = boQ^^^, containing g original monomers. View- 
ing the collapsing chain at a length scale smaller than 
R reveals straightened segments. However, the overall 
large scale structure still remains a random walk of seg- 
ments. This picture is clearly captured in snapshots from 
simulations, which we discuss below; see Fig. |^. 

If we neglect numerical factors associated with detailed 
shapes of the objects being pulled through the viscous 
liquid, the hydrodynamic equations of motion for coarse- 
grained monomers at scale R can be written: 



V - f/voAiR) 
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FIG. 2: Simulation snapshots of a collapsing A'^ = 512 chain. 
Time t is measured from the instant at which the solvent 
quality is changed from "good" to "poor". 



where rjo is solvent viscosity, / is the chain tension, and 
V is the drag velocity. The quantity A(R) measures the 
friction experienced by objects due to the solvent. For 
Zimm dynamics, A(R) depends on the largest geometric 
dimension of the object: A{R) ~ R. In Rouse dynamics, 
A{R) ~ g regardless of its geometric size. For now, we 
assume Zimm dynamics holds, and hence, the velocity 
scales as u ~ f /rjoR. 

We account for the chain tension / as follows. By 
bringing gi ^ L monomers forming a string of length 
L < R into a globule, one gains free energy AE/kBT — 
giAT/9. Thus, 



L hn 



(2) 



Note this force / is independent of scale R. In this esti- 



mation of force, we neglect the surface term correction, 
fci3T/&og-i/3(Ar/6')4/3, as it decays relatively rapidly 
with increasing g. There are also external forces act- 
ing on each segment from the remainder of the chain. 
But because the segments form a random walk, the vec- 
toral average of these external forces acting on the ends 
of any segment is zero. Hence, the initially given Gaus- 
sian structure remains through the coarse-graining, as 
the length of the segments increases. 

The velocity v of motion on each scale R is controlled 
by Eq. |l|. After a time t ~ R/v ~ i?^, in which all 
segments move a distance i?, the configuration on scale 
R is "relaxed." We can assume that all g monomers 
of a subchain collapse to a cylindrical configuration, or 
segment, of length R{t) and diameter d{t). The volume 
of the cylinder RdP should be approximately equal to the 
total volume of g monomers: V ~ glf^. Since the largest 
dimension R of the cylindrical object remains the same 
as for initial Gaussian chain R ~ 5^^^, we conclude that 
g t, and d ~ g^^^- Hence, we arrive at scalings R ^ t^^^ 
and d t^/^. (For Rouse dynamics, one obtains R ^ t^^^ 
and d ^ t^/^). 

It is important to realize that these cylindrical config- 
urations exist only in a statistical sense. Any uniformly 
dense "cylindrical" configuration would be immediately 
broken into an array of globules, connected by strings, 
due to the capillary instability, jist The relaxation time 
for the capillary effect t ^ d, including the coalescence of 
globules, is short, compared to streamlining time t ^ R^ 
(or t ~ d^). Thus, the formation of the pearl- necklace 
structure (within the cylinder) does not require modifi- 
cation of the streamlining description. Our cylinders are 
arrays of globules characterized by an average diameter of 
d and inter-globule distance d, and overall array length R. 
Mass conservation then requires that each globule con- 
sists of g'^^'^ monomers, and the total number of globules 
is 

Furthermore, there is no energy barrier associated 
in bringing two separated globules together. This is 
in contrast to models of other authors, who main- 
tain that globule interaction is the main mechanism of 
relaxation. p^, |l3j Instead of considering the processes 
of coalescence of globules, and the competition in growth 
and evaporation in detail, here we adopt the simple view- 
point that the process of streamlining on scale R results 
in the formation of statistical cylinders of size R. The 
contraction of each segment drags the globules to the 
larger-scale cylinder. The overall structure is a random 
walk of such cylinders, which still displays a Gaussian 
structure on larger scales until the end of coarse-graining 
stage. 

This mechanism allows direct derivation the time- 
dependent structure factor S{k, t) for a collapsing chain. 
S{k, t) is a useful function, because it can be both mea- 
sured experimentally and directly computed from simu- 
lations. Below we report the structure factor S{k) of a 
polymer chain made of cylinders, which in turn are made 
of linear necklaces of spheres. The structure factor can be 
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written as a product of structure factors of its elements: 



S{k,t) = Sgauss[Rit)]S,yl[Rit)]Sglb[dit)]. 



(3) 



S gauss {k, t) is the structure factor of Gaussian chain con- 
sisting of Ni{t) = N/g{t) unit cyHnder elements of aver- 
age length R{t). Therefore, 



S. 



gauss 



1 f 2/iVi [e-"(^i+i) - e-"] - + 



(1 - e-<^y 



(4) 

where a = \kR{t)] /3. For large fc, we note that this 
expression has a log-log slope of -2, as expected. Scyi 
and Sgib respectively are the structure factors of cylinder 
with length R{t) and a globule of size d{t): 



Scyl{R,k) = 

Sgib{d,k) = 



--— [cos kR + kRSi{kR) - 1] (5) 



3 1' 
— ^ (sin kd — kd cos kd) 



(kd) 



(6) 



Si(x) is the sine integral of x. The length of cylinders 
R{t) and the diameter of globules d{t) are the parameters 
which follow the coarse-graining rule suggested above. 
The large-fc log-log slopes of the relations in Eq. ^ and 
Eq. ^ are -1 and -4 respectively. 

As a test of our theory, we computed the structure fac- 
tors for a simulated polymer chain undergoing collapse. 
The collapse of a homopolymer chain in a poor solvent 
was simulated via molecular dynamics (MD). The pro- 
totype system is a single linear freely-jointed polymer 
chain of length A'' = 512 monomers, immersed in a sol- 
vent of 8x10"' simple particles, contained in a cubic box. 
The total particle number density, p, is 0.85 a~^, giv- 
ing a box dimension of 45.6 a. All particles interact via 
standard Lennard- Jones potentials. For "good" solvent 
conditions, the solvent-solvent, monomer- monomer, and 
solvent-monomer pair interaction, designated J7g(r,j ), is 
a shifted repulsive Lennard- Jones potential, cut off at 
r — 2^/^ a. For "poor" solvent conditions, the solvent- 
solvent and monomer-monomer pair interaction is the full 
Lennard- Jones potential, UL.j{rij), with a well depth e 
of 1.0 ksT, cut off at r = Tc = 2.5 a, while the solvent- 
monomer pair interaction is kept as the repulsive Ug . The 
shift s of Ulj is set such that Ul.i{tc) — 0. This sim- 
ple prescription was used previously in smaller-scale ho- 
mopolymer collapse simulations. [Q Particles connected 
to each other along the backbone of a chain are con- 
strained via a harmonic potential with an average length 
of 0.97 a and a standard deviation of 0.05 a. 

The equations of motion are integrated using the 
velocity-Verl et algorithm , with a time step At of 0.001 r, 
where r — ma'^ /ksT ^ and m = 1 . Energy is measured 
in units of /cbT, which prescribed that the particle di- 
ameter, do, be exactly a, for both monomers and solvent 
particles. Thus the total excluded volume fraction in the 
system is about 0.45. 

The procedure for a single collapse trajectory is as fol- 
lows. The initial configuration was generated by growing 



a random walk of length N = 512 and step length I = 
0.97 (7, biased such that 1-3 monomers automatically 
have no overlap. Then the remaining 80 000 solvent par- 
ticles are added at random locations in the box. Parti- 
cle overlaps are then removed using a previously detailed 
warmup technique, [p^ which incrementally increases the 
the excluded volume diameter of the particles from zero 
to its full value over about 5000 MD time steps in "good" 
solvent conditions. During this phase, a Langevin-type 
thermostat with friction F = 5.0t~^ is employed. The 
amount of time used for this warmup was kept as low as 
possible in order that the initial configurations remain as 
Gaussian as possible without significant swelling due to 
the actual presence good solvent conditions. 

Immediately following the warmup, designated for con- 
venience as t = 0, the quench into poor solvent conditions 
is performed via instantaneously switching the monomer- 
monomer and solvent-solvent pair interactions from Ug 
to Ulj- The MD integration is run for several thousand 
T, until the chain has fully collapsed into a single glob- 
ule. During this phase, no thermostat is employed, and 
the temperature was observed for all runs to fluctuate 
around 1.0 by less than 0.1%. The collapse is monitored 
using the radius of gyration of the chain as a function 
of time. During this integration, configuration snapshots 
are stored every ten thousand time steps (10 r) for later 
analysis. All simulations were run in parallel on 256 pro- 
cessors using domain decomposition. [0 

A total of 40 independent collapses were simulated, 
each beginning with a unique chain conformation. The 
time-dependent structure factor, S{k,t), was computed 
at times t = 0, 10, 20, 40, 80, and 160 r, by averaging 
the static structure factor over all chains at each of those 
times. We note that the standard deviations were about 
20% of the average, which is consistent with large vari- 
ances in size among the different initial Gaussian chains. 

The primary result of this work is a comparison of 
S{k,t) predicted by Eq. |^ and S{k^t) computed from 
the MD simulations. In order to make this comparison 
quantitative, we introduced the necessary prefactors in 
the scaling relations R ~ t^/"^ and d ~ t^/'^. We then 
performed a least-squares, conjugate gradient minimiza- 
tion ||l^ on the sum of squared errors between the theo- 
retical and computed S{k,t) for all points for fc < 1 and 
for all times simultaneously. The resulting comparison is 
shown in Fig. ^. Note that the theoretical and computed 
S{k, t) agree remarkably well. In particular, both display 
the three asymptotic slopes of the structure factor com- 
ponents discussed previously. The theory and simulation 
are linked through the following relations: 

R[t) = At^/^ (7) 
d{t) = Bt^/^ (8) 

with A = 1.88(9) cr/Ti/2 and B = 0.78(8) g/t^/^. The 
poorest agreement appears for the latest time, perhaps 
signifying the end of the coarse-graining stage of the col- 
lapse. Due to the large inherent variability of the simu- 
lation data, we do not report goodness-of-fit statistics. 
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FIG. 3: The single ehain time-dependent structure factor, 
S{k,t). Comparison between theoretical predictions (lines) 
and simulation (points). The t = structure factor is mea- 
sured over all chains immediately prior to the quench, and 
is best represented as S ~ k~^'^, close to the theoretical 

To complete our picture, we discuss the late-stage 
collapse dynamics. After the coarse graining stage, a 
chain becomes a single "cylinder", where the globules 
are aligned linearly. The terminal globules experience 
net forces due to the finite size of the chain. The late 
stage is therefore characterized by "end effects." The 
typical late stage consists of two globules with connect- 
ing string between these two. In this regime, the dynam- 
ics are controlled by tension between the two globules 
and the hydrodynamic friction experienced by the pearls 



(not by the strings). Both terminal globules move to- 
ward the common center of mass by "eating" the string. 
This procedure is deterministic. Considering the hydro- 
dynamic friction for the sphere with size Rp is ^ GirrfRp 
and Rp ~ N^^^, the typical distance between the termi- 
nal globules should be N^^'^ (reminiscent of a Gaussian 
fractal) , one can estimate the average characteristic time 
for the late stage as t; - R/V ~ N^^^ {rf^ ~ N^^^ for 
Rouse dynamics). Future work will offer a quantitative 
comparison of this picture with that observed in ongoing 
MD simulations. 

Wc have presented a ncnv tlic;ory c;apable of describing 
the dynamics of polymer collapse in poor solvent. We 
have proposed that the initially fractal chain recursively 
straightens on length scales R{t), forming an increas- 
ingly coarse-grained random walk of statistical cylinders 
of length R ~ t^/^ and diameter d ~ f^/^, in accordance 
with hydrodynamic equations of motion. These cylin- 
ders are statistical objects, whose underlying structure 
consists of necklace-like configurations of randomly po- 
sitioned globules connected by strings. This theory suc- 
cessfully incorporates the capillary instability into a hy- 
drodynamic framework. The theoretical time-dependent 
structure factor is strongly supported by those computed 
from large-scale MD simulations. 
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